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Abstract. An increasing number of model results suggests that chiral symmetry is broken inhomogeneously 
in a certain window at intermediate densities in the QCD phase diagram. This could have significant effects 
on the properties of compact stars, possibly leading to new astrophysical signatures. In this contribution 
we discuss this idea by reviewing recent results on inhomogeneous chiral symmetry breaking under an 
astrophysics-oriented perspective. After introducing two commonly studied spatial modulations of the 
chiral condensate, the chiral density wave and the real kink crystal, we focus on their properties and their 
effect on the equation of state of quark matter. We also describe how these crystalline phases are affected by 
different elements which are required for a realistic description of a compact star, such as charge neutrality, 
the presence of magnetic fields, vector interactions and the interplay with color superconductivity. Finally, 
we discuss possible signatures of inhomogeneous chiral symmetry breaking in the core of compact stars, 
considering the cases of mass-radius relations and neutrino emissivity explicitly. 

PACS. 26.60.-c Nuclear matter aspects of neutron stars 


1 Introduction: inhomogeneous phases in 
dense quark matter 

Compact stars provide the only known realization of ultra- 
dense matter in nature. While the external layers of these 
stars are expected to be formed by more conventional 
forms of matter, like nuclei, neutrons, protons, and elec¬ 
trons, the composition of matter in the center and its 
properties are still unresolved. This in turn implies that 
observations of compact stellar objects might help deter¬ 
mine the properties of strong interactions at finite den¬ 
sities, which is one of the most challenging tasks of con¬ 
temporary nuclear physics. While ab-initio lattice calcula¬ 
tions are still unable to provide definite answers, the phase 
structure in this region might be extremely rich. Already 
half a century ago, it was speculated that the central den¬ 
sities of compact stars could be high enough to liberate 
quark degrees of freedom m- More recently this idea 
was challenged by the discovery of two compact objects of 
two solar masses m, requiring a rather stiff equation of 
state. This puts severe constraints on the appearance of 
quark matter in compact stars, but does not rule out this 
possibility [5] (see [5] in this topical issue for a detailed 
analysis and [7] for an overview). 

In this contribution we assume that quark matter is 
present in compact stars and discuss its properties and 
possible phenomenological consequences. In particular, in 
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the past decade several studies have suggested that dense 
quark matter might form crystalline, i.e. spatially inho¬ 
mogeneous structures, leading to a significant revision of 
the current picture of the phase diagram of quantum chro¬ 
modynamics. These structures are related to the conden¬ 
sation of fermions into pairs carrying a finite momentum, 
translating into an explicit spatial dependence of such con¬ 
densates. 

Prominent examples of inhomogeneous phases in dense 
quark matter are related to the phenomena of color super¬ 
conductivity and chiral symmetry breaking. Inhomoge¬ 
neous color superconductors might appear in isospin im¬ 
balanced quark matter, where the separation between the 
quark Fermi surfaces makes the formation of diquark pairs 
with finite momentum the energetically favored channel. 
Inhomogeneous chiral condensates on the other hand arise 
in the region between the low-density phase where chiral 
symmetry is spontaneously broken and the high-density 
restored one. In the usual picture limited to homogeneous 
chiral condensates made of quark-antiquark pairs, the pres¬ 
ence of a chemical potential induces a stress on the for¬ 
mation of these pairs, eventually leading to chiral restora¬ 
tion. On the other hand, at finite densities the forma¬ 
tion of pairs made of quarks and holes at the Fermi sur¬ 
face becomes a competitive condensation channel. Since 
pairs with large relative momenta are disfavored, the pre¬ 
ferred pairing mechanism leads to a condensate with a 
non-vanishing net total momentum, i.e. a spatially inho¬ 
mogeneous one (for a more detailed discussion, see e.g. 
0 ). 
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The theoretical investigation of inhomogeneous con¬ 
densates is mostly performed through effective models, 
most notably the Nambu-Jona-Lasinio (NJL) and Quark- 
Meson models. Most studies of this kind find an inhomo¬ 
geneous “island” at low temperatures and finite densities, 
like in the phase diagram depicted in fig. crystalline chi¬ 
ral condensates are thermodynamically favored in a finite 
chemical potential window, typically covering the usual 
first-order phase transition which is expected to occur for 
homogeneous matter. More recently, this picture was con¬ 
firmed by a Dyson-Schwinger calculation, directly based 
on QCD [9]. The actual extension of the inhomogeneous 
window is still uncertain: the Dyson-Schwinger and some 
NJL studies suggest that, at least at low temperatures, 
inhomogeneous chiral symmetry breaking might occur in 
a very wide range of densities wm, while other models 
predict a smaller size. In this sense, fig. could be seen as 
a rather conservative estimate. 

At any rate, the realization of inhomogeneous chiral 
symmetry breaking appears to be a rather robust predic¬ 
tion for the region of the phase diagram which is expected 
to be relevant for the description of compact stars. If quark 
matter is present in the core of compact stars, the forma¬ 
tion of crystalline condensates might therefore have sig¬ 
nificant consequences on the physical properties of these 
objects, possibly leading to new signatures. Since possible 
observables related to crystalline color superconductivity 
have already been discussed in m , here we will mainly fo¬ 
cus on inhomogeneous chiral condensates. Also, since our 
main objective is to discuss the relevance of inhomoge¬ 
neous quark matter for compact stars properties, we will 
mostly consider vanishing temperatures and try to present 
established results under a new, more astrophysics-oriented 
perspective. The interested reader can find a complemen¬ 
tary review with a more general discussion on inhomoge¬ 
neous chiral condensates, including their finite-temperature 
properties and resulting phase structures in [12] (see also 
m for a brief historical review). 

In the following, after describing basic features of in¬ 
homogeneous chiral symmetry breaking and its effects on 
some thermodynamical quantities of a dense quark sys¬ 
tem, we will discuss how the physical conditions realized 
inside compact stars, such as charge neutrality or the pres¬ 
ence of strong magnetic fields, might affect the properties 
of these crystalline condensates. 


2 Properties of inhomogeneous chiral 
symmetry breaking 

In this section we want to give a brief overview about the 
basic features of inhomogeneous chiral condensates and 
how they are formally derived in model calculations. To 
be specific, we consider an NJL model, defined by the 
Lagrangian 

>Cnjl = V' - m) Ip+ G {{ip'tpf' + , 

( 1 ) 



p (MeV) 

Fig. 1. Phase diagram in the /r-T plane. The shaded region in¬ 
dicates the inhomogenous phase, covering the first-order phase 
boundary (blue solid line) between the homogeneous chirally 
broken and restored phases. This phase diagram is a reproduc¬ 
tion of the result shown in ref. [14] and was obtained within 
an NJL model with the RKC ansatz, see sect.[^ 


where ip denotes a quark field with bare mass m, inter¬ 
acting by local four-point vertices proportional to the di¬ 
mensionful coupling constant G. For simplicity, we be¬ 
gin by considering here the standard two-flavor model 
with scalar-isoscalar and pseudoscalar-isovector interac¬ 
tions |15] . while we will discuss various model extensions 
in later sections. 

In order to analyze the thermodynamic properties at 
temperature T and quark chemical potential /i = 

(fiB being the baryon chemical potential) we need to cal¬ 
culate the grand potential per unit volume V, 

I2(T,p) = -^logZ(r,p), (2) 

where 2 denotes the grand canonical partition function. 
To this end, we perform a mean-field approximation in the 
presence of the scalar and pseudoscalar condensates 

S(x) = (ipip) and F’(x) = (ipi'y^T^ip), (3) 


respectively. Here we assume that only the third isospin 
component of the pseudoscalar condensate is non-vanish¬ 
ing. By this, we exclude the possibility of charged pion 
condensation which at any rate is not expected to occur in 
compact stars, where the isospin chemical potential stays 
well below the pion mass (sect. 3.1). 

We also assume that the condensates are time indepen¬ 
dent. However, since we are interested in inhomogeneous 
phases, we must retain their dependence on the spatial co¬ 
ordinate X. As a consequence, the derivation of the mean- 
field thermodynamic potential becomes significantly more 
difficult than for space-independent condensates. As dis¬ 
cussed in detail in ref. [El, one finds 


1^Mf(F, /r; H) — Ufcm T G^ond : (d) 

with the condensate part 

G,ond = ^ J^d^xG{S^{x)+P^{x)) (5) 
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and the kinetic part 


The mass function reads 




1 -b e T 


( 6 ) 


where A labels the eigenstates of the effective Dirac Hamil¬ 
tonian 


H(x) = 7 ° [-ifd, + m- 2G(5(x) + i-f\^P{x))] , (7) 

and E\ are the corresponding eigenvalues. In contrast to 
the homogeneous case, the diagonalization of H for arbi¬ 
trary space-dependent condensate functions is rather dif¬ 
ficult and must in general be performed numerically. After 
that, rnust be minimized with respect to the func¬ 

tions 5'(x) and P(x) in order to hnd the energetically fa¬ 
vored ground state. Since this problem has not yet been 
solved in 3 -b 1 dimensions, one considers simple ansatze 
for the condensate functions. This reduces the problem to 
finding the minimum of 17 mf with respect to a limited 
set of variational parameters. Most parametrizations used 
in literature correspond to one-dimensional spatial mod¬ 
ulations. Two- and higher-dimensional modulations have 
been studied as well, but (unlike for color superconduc¬ 
tors) were found to be disfavored against one-dimensional 
ones, at least within the mean-field approximation nani]. 

In the following we will restrict ourselves to the chiral 
limit (to = 0) and consider two one-dimensional ansatze 
which even allow for an analytical diagonalization of H 
and are therefore the most commonly used in the litera¬ 
ture. Choosing the modulation to be along the ^-direction 
and introducing the complex mass function 

M{z) = -2G{Siz)+iP{z)), ( 8 ) 

the most simple (and therefore most popular) one is the 
so-called (dual) chiral density wave (CDW) [THIIT^ 


M(z)=Z\e*«^ (9) 

i.e. a single plane wave with amplitude A and wave num¬ 
ber q, corresponding to a spatial period of L = 2Tr/q. 
When moving along the ^-direction, the CDW rotates 
around the chiral circle, transforming scalar condensates 
into pseudoscalar ones and vice versa: 

Siz) = -^ cos (qz), P{z) =-^ sin (qz). (10) 

A more sophisticated ansatz is the “real kink crys¬ 
tal” (RKC), corresponding to an array of domain-wall 
solitons. It was derived from the analytically known so¬ 
lutions of the 1 -b 1 dimensional Gross-Neveu model [20] . 
embedding them into 3 -b 1 dimensions by making use 
of Lorentz boosts into the transverse directions [Hj. Al¬ 
though there is no general proof, the RKC modulations 
have so far turned out to be the energetically most fa¬ 
vored ones among all ansatze considered in the literature, 
including higher-dimensional ones uniiiz], so that to our 
current knowledge they consitute the most likely shape 
realized in the inhomogeneous island. 


M{z) = A^sn{Az\v) (11) 


and depends on two parameters, A and iz. Here sn(a|iz) 
denotes a Jacobi elliptic function with elliptic modulus v G 
[0,1], interpolating between tanh(a) for zz = 1 and sin(a) 
for zz = 0. The amplitude of the modulation is then given 
by A^/l/ and its spatial period by L = 4K(zz)/Z\, where 
K(zz) is a complete elliptic integral of the first kind. In 
contrast to the CDW, the RKC is real (hence the name), 
i.e. the pseudoscalar condensate vanishes identically. 

Examples for the shape of this mass modulation are 
displayed in fig. In the lower part of the hgure we also 
show the corresponding density profiles, which have been 
derived in ref. HI]. As one can see, the structure is most 
pronounced near the phase boundary to the homogeneous 
phase with broken chiral symmetry. Here the density is 
sharply peaked around the zero crossings of the mass func¬ 
tion, while it is suppressed in the regions where |M(z)| is 
large. Increasing the chemical potential, the mass ampli¬ 
tude decreases (see also hg. and the density profiles 
become more and more washra out, smoothly approach¬ 
ing the corresponding density in the restored phase. 

This observed behavior suggests us that the shape of 
the density is largely determined by the mass function. 
In order to quantify this, we introduce a local Fermi gas 
(LFG) approximation, where the density at each point in 
space is calculated by the corresponding expression for a 
homogeneous ideal Fermi gas of quarks with mass M{z), 
i.e. at T = 0, 

ULFoiz) = ^ - \Miz)\^f^ 9i^,-\Miz)\). (12) 


The resulting densities are shown in the lower panel of 
fig- i as dotted lines. Comparing them with the exact 
mean-field results from HI] (solid lines) we can see that, 
although the two share a qualitatively similar behavior, 
the actual density for the RKC solutions is always larger 
than the LFG one. In particular, the RKC density over¬ 
shoots the density for the restored phase HUES], while 
riLFG/nrest Is always < 1, as can easily be seen from 
eq. @. 

While for the RKC the density is spatially varying, 
following the changes in the mass function, in the case 
of a CDW modulation, for which \M{z)\ = A = const., 
the density is constant in space. Nevertheless, similar to 
the RKC case, the density for the CDW is larger than the 
one obtained for a spatially homogeneous order parameter 
M = A. 


In the numerical examples presented in this section, 
we employed a Pauli-Villars type regularization (for more 
details, see M), and the model parameters are fitted to 
reproduce the chiral-limit value = 88 MeV of the pion 
decay constant and a constant constituent quark mass 
of 300 MeV in vacuum. Another popular regularization 
scheme, which is often used when dealing with magnetic 
fields, is the Schwinger proper time method, for which 
the parameters are usually fitted to a slightly larger con¬ 
stituent mass of 330 MeV, see e.g. [I^[23] . We note in any 
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Fig. 2. Mass functions (upper row) and density profiles normalized to the density of chirally restored matter (lower row) for an 
RKC modulation, evaluated at zero temperature and = 308 (left), 325 (center) and 345 MeV (right), corresponding to regions 
near the phase boundary to the homogeneous chirally broken phase, in the center of the inhomogeneous phase and close to the 
phase boundary to the restored phase, respectively. For the density, the actual RKC results (solid lines) are compared with the 
ones obtained in the LFG approximation (dotted lines). 


case that the choice of a different regularization does not 
qualitatively alter the results (with the notable exception 
of a sharp three-momentum cutoff, which is not appropri¬ 
ate when dealing with inhomogeneous condensates since 
quark momenta are not conserved quantities). 

Minimizing the thermodynamic potential for the RKC 
ansatz with respect to A and v at given values of ^ and T 
then gives rise to the phase diagram shown in fig.j^ which 
we already discussed in the Introduction. Using the CDW 
ansatz and minimizing with respect to A and q yields 
a rather similar phase diagram. Although the CDW is 
disfavored against the RKC [ll] , working with the simpler 
CDW ansatz is therefore often sufficient to get a rough 
picture of the phase structure for inhomogeneous phases. 

The main qualitative difference between the two modu¬ 
lations concerns the transition from the homogeneous bro¬ 
ken to the inhomogeneous phase, which is first order for 
the CDW and second order for the RKC. This can also 
be seen in fig. where the amplitudes and periods are 
displayed for both modulations. At the phase transition, 
the amplitude drops discontinuously for the CDW ansatz, 
while it is continuous for the RKC. For the latter the pe¬ 
riod L diverges at this point, enabling a smooth transition 
to the homogeneous phase. The phase transition from the 
inhomogeneous to the restored phase, on the other hand, is 
second order in both cases, with the amplitudes smoothly 
going to zero. 

Except for the region close to the phase transition to 
the homogeneous broken phase, the periods are almost 
identical for both ansatze, roughly being of the order of 
2-3 fm. The amplitudes, on the other hand, differ by 
about a factor of 2. These results are consistent with pre¬ 
dictions of a Ginzburg-Landau analysis near the second- 
order transition to the restored phase 




p (MeV) 


Fig. 3. Amplitudes (upper panel) and spatial periods (lower 
panel) of the mass functions M{z), allowing for CDW (dashed) 
or RKC (dash-dotted) modulations, as functions of the quark 
chemical potential at T = 0. 


The spatially averaged quark-number densities 
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Fig. 4. Space-averaged quark-number density h as a function 
of the quark chemical potential at T = 0. The dashed and 
dash-dotted lines indicate the results obtained when allowing 
for CDW and RKC modulations, respectively. The density for 
the restored solution (M = 0) is indicated by the solid line. 
Note that the baryon number densities ub are related to the 
quark ones by ub = njj,. 

are displayed in fig. together with the corresponding 
density for the restored solution. At fixed chemical poten¬ 
tial, h is smallest for the RKC and largest for the restored 
solution, with the CDW lying in between. Again, this can 
qualitatively be explained as a suppression effect due to 
the effective quark mass, which vanishes in the restored 
phase and is largest in the RKC[^ 

For the RKC ansatz, the inhomogeneous phase starts 
essentially at ft = 0, corresponding to infinitely sepa¬ 
rated domain-wall solitons. For the CDW solution, on the 
other hand, the minimum density is around 1.4 no, where 
no « 0.5 fm“^ is the quark number density in symmetric 
nuclear matter at saturation. In both cases we find the 
upper density of the inhomogeneous island to be about 
2.2 no. Keeping in mind that these numbers are model and 
parameter dependent and could also be higher, this lies 
in the ballpark of densities expected in a possible quark 
phase in compact stars. 

The pressure p = —(17 — f2vacuum) of the various so¬ 
lutions at T = 0 as function of p is shown in fig. As 
usual, the favored solution is the one with the largest pres¬ 
sure. At low chemical potentials, this is the homogeneous 
broken solution (dotted line). Excluding the possibility of 
inhomogeneous phases, there would be a first-order phase 
transition to the restored phase (solid line), characterized 
by the crossing of the two lines at the critical chemical 
potential. However, by allowing for inhomogeneous solu¬ 
tions, we see that in this regime the RKC (dash-dotted 
line), which sets on in a second-order phase transition at 
a slightly lower chemical potential, is more favored. Ad¬ 
ditionally, we observe that the CDW is always disfavored 
against the RKC, but still favored over the homogeneous 
solutions in most of this region. At higher chemical poten¬ 
tials the pressure curves of the inhomogeneous solutions 
approach the one of the restored phase and eventually 

^ This remains true when we take the spatial average of 
\M{z)\\ 



Fig. 5. Pressure at T = 0 as a function of the quark chem¬ 
ical potential for different mass functions. Solid: restored so¬ 
lution, dotted: homogeneous broken solution, dashed: CDW, 
dash-dotted: RKC. 

reach it in a second-order transition (not shown in the 
figure). 

Finally, in fig. we compare the corresponding equa¬ 
tions of state at T = 0, i.e. the pressure as a function of 
the spatially averaged energy density, 

e = —p -\- pn. (14) 

At the first-order phase transitions from the homogeneous 
broken to the restored phase or the CDW, h (and hence 
e) are discontinuous, leading to horizontal jumps in the 
equation of state. The transition to the RKC, on the other 
hand, is continuous. However, since the density grows very 
quickly above the transition point (cf. fig.|^, the resulting 
equation of state is very similar to the discontinuous ones. 
In general, the inhomogeneous equations of state and in 
particular the RKC one are stiffer than the homogeneous 
one, but the differences are rather small. We therefore 
anticipate that the effect of inhomogeneous phases on the 
mass-radius relation of compact stars will be negligible. 

3 Including further elements of compact-star 
physics 

After discussing the basic properties of inhomogeneous 
chiral condensates within a rather simple model, we now 
study the effects on them of additional physical conditions 
encountered in compact stars. 

3.1 Electric neutrality 

In the previous section we have assumed that up and down 
quarks share the same chemical potential p. In compact 
stars, however, one has to take into account weak decays 
and ensure global electric neutrality. In order to describe 
this situation, our model must be extended to contain lep¬ 
tons. Since those are only weakly coupled to the quarks, 
the total thermodynamic potential can be decomposed 
into a quark and a leptonic part, 

^tot = 


(15) 
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Fig. 6. Equation of state (pressure p vs. spatially averaged en¬ 
ergy density e) at T = 0 for homogeneous matter (solid) and in 
the inhomogeneous phase with CDW (dashed) or RKC (dash- 
dotted) modulations. All equations of state begin at p = e = 0 
in the homogeneously broken phase. The dotted horizontal 
lines indicate discontinuities related to first-order transitions 
to the restored phase or to the CDW solution. The inset shows 
an enlarged detail of the transition region. 


Moreover, neutrinos can freely leave the star, so that we 
only need to take into account charged leptons. 

This system is characterized by two conserved global 
charges: the net quark number and the total electric charge. 
Therefore f2tot depends on two independent chemical po¬ 
tentials, the quark number chemical potential ^ and the 
electric charge chemical potential /rg. The condition for 
(spatially averaged) electric neutrality is then given by 


_ df2tot ! 


( 16 ) 


The chemical potentials of the individual particle speci¬ 
es are given as linear combinations of p and /ig , according 
to their quantum numbers, i.e. 


2 1 , , 

Ai« = Ai+3MQi ^J-d = ^J■s= ^^J-Q (17) 


for the quarks and 


= -MQ 


(18) 


for the (negatively charged) leptons. 

Neglecting Coulomb effects aside from the neutrality 
constraint, we can approximate the leptonic part of the 
thermodynamic potential by an ideal gas. The remaining 
problem is then to calculate the quark part at nonzero 


/ig and to determine /ig from eq. (161. 

In the quark sector, if we restrict ourselves again to up 
and down quarks, /ig is equivalent to an isospin chemical 
potential /i/ = /i„—/id = /ig. Inhomogeneous phases in the 
presence of a nonvanishing isospin chemical potential have 
been discussed in refs. |25LI26j within a Ginzburg-Landau 
analysis, which is expected to be valid near the chiral 
critical point but not reliable at the low temperatures of 
compact stars. Recently, inhomogeneous phases in isospin 
asymmetric matter have also been studied within the NJL 
model [27]. Starting from the Lagrangian of eq. (Q a 


mean-field approximation was performed, considering the 
flavor-diagonal scalar and pseudoscalar condensates 

'5'/(x) = (//) and P/(x) = (/nV), (19) 

/ € {m, d}, and a generalized CDW ansatz of the form 


Sf{z) 


AG 


cos (qfz ), 


Pfi^) 


AG 


sin (qfz). 


( 20 ) 


Note that the condensates defined in eq. @ are related 
to the above ones hy S = Su + Sd and P = Pu — Pd, so 
that eq. (101 is recovered if we choose = Ad = A and 
qu = -Qd = q- 


For qu = —qd and arbitrary amplitudes Au and Ad 
the effective Dirac Hamiltonian of this ansatz can still be 
diagonalized analytically. Minimizing the thermodynamic 
potential for this case it was found that the inhomoge¬ 
neous phase shrinks with increasing |/^/|, but it still sur¬ 
vives at values of /r/ needed to obtain electrically neutral 
matter [27] 

It is plausible, however, that the assumption = —qd 
is too restrictive: As we have seen earlier, the spatial pe¬ 
riodicity and, hence, the wave number is a function of the 
chemical potential. At nonzero /r/, where the chemical po¬ 
tentials of up and down quarks differ, we should therefore 
expect that unequal values of |(/ii| and |( 7 d| could be fa¬ 
vored. Unfortunately, for this case the eigenvalues of the 
Hamiltonian are not known analytically, so that H has to 
be diagonalized numerically. Obviously, this complicates 
the calculations considerably. In particular, for the nu¬ 
merical diagonalization a fixed ratio R — qu/qd of the 
wave numbers must be considered in order to be able to 
map the momenta onto a discrete grid. While in princi¬ 
ple any real number can be approximated by a given R 
and therefore all possible values for and qd can be ex¬ 
plored, in practice of course only a few selected ratios can 
be considered. 

In ref. [IT] the analysis was therefore restricted to the 
values R G { —1,|,2,5}, and it was determined which 
of these ratios is the most favored one at fixed |/r/| = 
60 MeV. The resulting phase diagram is shown in fig. [^ 
Here the phase structure is shown in the plane spanned by 
fd = + Md) = Id + \ldQ and the temperature T. One 

can see that different ratios are favored in different areas 
of the inhomogeneous phase, with the general tendency 
that R decreases with increasing p,. In a more realistic 
scenario one would of course expect that R varies con¬ 
tinuously with p and T. However, even within this very 
limited analysis, we find that the size of the inhomoge¬ 
neous region at nonvanishing /i/ is roughly the same as 
for fii = 0 (cf. fig. [^, so that this region might become 
even bigger if more ratios are considered. 

So far we restricted ourselves to the standard NJL- 
model Lagrangian, defined in eq. 0- However, for asym¬ 
metric quark matter, the exact isospin structure of the 


^ Solving eq. (161 for this ansatz, /tj was found to be around 
70 — 80 MeV in the inhomogeneous phase at T = 0. This is 
roughly the same result one gets for the restored solution in 
this regime. 
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Fig. 7. Phase diagram in the plane spanned by the averaged 
chemical potential fi = |(Mii + fJ-d) and temperature T at fixed 
isospin chemical potential |/ii| = 60 MeV for the generalized 
CDW ansatz eq. (20l with R — qu/qd £ { —1,|,2,5}. The 
shaded area indicates the region where an inhomogeneous so¬ 
lution is favored over homogeneous ones. Different colors cor¬ 
respond to different values of the most favored R, as indicated 
by the labels. The figure was taken from ref. [ 13 - 


interaction becomes relevant. In particular, eq. Q can be 
viewed as a special case of the more general Lagrangian [28l 


C= 

( 21 ) 


with the current quark matrix m = diag^(m„, m^;) and 
two SU (2) L X SU (2)ij X C/y (1) symmetric interaction terms, 
proportional to the coupling constants Gi and G 2 , respec¬ 
tively. The first term respects an additional UaI^) sym¬ 
metry, which is explicitly broken by the second. 

Considering again the condensates defined in eq. (19), 
one finds that the mean-field Hamiltonian depends on the 
effective mass operators 


Mf=mf-AGi {Sf + Pf - 4G2 (Sh - Ph ( 22 ) 

with the flavor indices f,h G {u,d}, f ^ h. From this 
expression we can see that the Gi term of the interac¬ 
tion dresses quarks of flavor / with condensates of the 
same flavor, while G 2 leads to mixing. In particular in the 
standard NJL Lagrangian eq. ([^, which is recovered for 
Gi = G 2 = G/2, the condensates of both flavors con¬ 
tribute equally to both masses, so that for rriu = md, the 
effective mass operators and Md are always equal as 
well. In this sense eq. ([^ corresponds to the maximally 
flavor-mixed case. 

For G 2 = 0, on the other hand, the flavors decouple. 
In this case the quark part of the thermodynamic poten¬ 
tial becomes Gq = -|- f2d, where Gf = only 


depends on the chemical potential The equation of 
state at nonzero fj,Q as a function of ^ is then trivially 
obtained from the one at /xq = 0 through eq. (171. In par¬ 
ticular the phase boundaries related to the properties of 
the up and down quarks are shifted relative to each other 
by fiQ, so that there are now regions where one flavor is 
homogeneous while the other one is inhomogeneous. 


In ref. [29] a realistic value for the ratio G 2 /G 1 was 
estimated to be around 0.2, i.e. the flavors are neither 
decoupled nor maximally mixed. A detailed discussion of 
this case will be given in ref. [30] . 


Finally we note that instead of considering a homo¬ 
geneous leptonic background and imposing the neutral¬ 
ity condition for the spatially averaged charge density, we 
could in principle do better and allow for an inhomoge¬ 
neous lepton density as well. While the density gradient 
would increase the kinetic energy of the electrons, it could 
lower the local Coulomb energy, and the optimal shape 
would result from balancing these two contributions. This 
possibility has not yet been considered. However, given the 
large incompressibility of electrons, a considerable varia¬ 
tion of their density over the length scale of the quark 
structures is unlikely, so that the assumption of a con¬ 
stant electron background is probably not too bad0 


3.2 Magnetic fields 

A very important element characterizing compact stars 
is the presence of strong magnetic fields. For magnetars, 
surface fields up to lO^^G have been measured, and even 
higher values are expected to be reached in their cores, 
possibly reaching up to 10^®G m- More recently even 
higher values have been suggested |32], although under 
such extreme regimes large deformations are expected to 
be present and the star structure is likely destabilized (33] 

[Ml- 

Aside from the well known effects on the usual ho¬ 
mogeneous chiral symmetry breaking mechanism, such as 
the phenomenon of magnetic catalysis |331I331[37] (see [351 
139] for recent reviews), the presence of strong magnetic 
fields has a significant influence on the formation and the 
properties of inhomogeneous phases. To see this, we start 
once arain from the simple two-flavor NJL Lagrangian 
of eq. (IH, and include the coupling with a background 
magnetic field by the replacement 9^ —)■ + 

iQAf^, where is the electromagnetic field and Q = 
diag(e„,ed) = diag(|e, —|e) the electric charge matrix 
in flavor space (e being the unit electric charge). Due to 
their unequal charges, the interaction with the magnetic 
held distinguishes between quark havors, effectively break¬ 
ing isospin symmetry. 

For a GDW-type ansatz (eq. ([^) in the presence of a 
static background magnetic held H pointing in the z di¬ 
rection, the quark energies for each havor have been cal- 


® We thank the referee for pointing this out to us. 
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culated in [23] and are given by 

+Pz+ g) +2|e/i7|n n>0. 

(23) 

By inspecting this expression one can see that the quark 
momenta transverse with respect to the direction of JI 
become quantized into ^ 2n|e/i?|, with the integer n 
labelling the so-called Landau levels. In turn, this implies 
that at the lowest Landau level (LLL), for which n = 
0, the problem becomes effectively l-|-l-dimensional|^ a 
configuration in which inhomogeneous chiral condensation 
is always favored gnEi]. 

In particular as a consequence of the spectral asymme¬ 
try at the LLLjj at low temperatures the CDW becomes 
the thermodynamically favored state over homogeneous 
matter for any p > 0. As can be seen from fig. [^ at low p 
the system is characterized by a small but nonzero value of 
the wave number q, until a phase transition is reached in 
the chemical potential region where the inhomogeneous 
phase would appear even for H = 0. In the low p re¬ 
gion the value of q is strongly influenced by the magnetic 
field, mainly due to the presence of an anomalous term 
^ \eH\fiq in the thermodynamic potential. Therefore the 
effect is more visible in the lower panel of fig. [^ corre¬ 
sponding to a larger value of H, although it also exist 
for smaller magnetic fields (upper panel). After the phase 
transition, the value of q becomes comparable with the 
one obtained in absence of a magnetic field, suggesting 
that the dominant physical mechanism there is the usual 
particle-hole pairing at the Fermi surface. 

If the magnetic field is sufficiently large, the amplitude 
of the CDW does not melt completely as p increases and 
the inhomogeneous phase extends to higher densities (for 
a more detailed discussion on the behavior of the CDW 
order parameters in a magnetic-field background we re¬ 
fer the reader to [13]). This suggests that chiral restora¬ 
tion might never be reached, although the enhancement 
of inhomogeneous phases at high densities might at least 
partly be a model artifact m- At high densities we ex¬ 
pect in any case color superconductivity to be the favored 
state for quark matter, see sect. |3.4] 

Being energetically favored, the CDW phase is char¬ 
acterized by a higher pressure when compared to the ho¬ 
mogeneous broken and restored ones, as can be seen from 
fig-0 We note at any rate that this effect is significant 
only for relatively large values of the magnetic field, while 
it becomes almost negligible for H < lO^^G. 

While so far we found that in 3+1 dimensions the RKC 
is favored over the CDW (cf. fig. [^, the introduction of 
a magnetic field and the related dimensional reduction at 
the LLL might change this picture. In particular, since in 
1+1 dimensional systems complex modulations are typ¬ 
ically favored over real ones mi, it is possible that the 

^ Indeed, the n = 0 eigenvalues are exactly the energies of 
a quark in a CDW-type condensate background in the chiral 
Gross-Neveu model. 

® For a more detailed discussion, see [42] . 




p (MeV) 


Fig. 8. Order parameters in a magnetic field for a CDW 
modulation as a function of the quark chemical potential, at 
given magnetic fields intensities. Top: eH = (28 MeV)^ {H = 
4.4 X 10^'^G), bottom: eH = (90 MeV)^ {H = 4.5 x 10^®G). 
Similar results have been presented in ref. |23| . 

CDW reappears in the phase diagram when a magnetic 
field is present. Moreover, it has been suggested in [33] 
that for finite H a hybrid chiral condensate (HCC) akin 
to the twisted kink crystal introduced in mi might be¬ 
come thermodynamically favored. This particular modu¬ 
lation can be interpreted as a generalization of the RKC 
to include a phase factor, or equivalently as a generaliza¬ 
tion of the CDW to include a spatially varying amplitude, 
and is given by 

Mhcc{z) = Z\v^sn(Z\z|:^)e®*^ . (24) 

For (7 = 0 this solution obviously reduces to the RKC, 
and if (/ = 1 it becomes thermodynamically degenerate 
with the CDW. If both g —>■ 0 and z/ —?> I, it approaches 
the homogeneous limit. While in light of the results dis¬ 
cussed in sect. [^ we would expect the imaginary part of 
this modulation to be zero, now the magnetic field always 
favors nonzero values for q, so that the order parameter 
is always complex. This implies that in the region of the 
phase diagram where the RKC inhomogeneous island used 
to appear, the favored solution is now given by the HCC, 
while in the other regions this ansatz reduces to the CDW 
modulation, except when chiral symmetry is restored. The 
p — H phase diagram at zero temperature looks then like 
fig. |10| for nonzero magnetic fields the low-density region 
is occupied by the CDW phase with small q, which then 
turns into the HCC after a critical chemical potential is 
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Fig. 9. Pressure at T = 0 in a background magnetic field as 
a function of the quark chemical potential for homogeneous 
matter and a CDW modulation. Top: eH = (90 MeV)^ {H = 
4.5 X 10^®G) bottom: eH = (280 MeV)^ {H = 4.4 x 10^®G) 

reached. The region of large /r and low H is chirally re¬ 
stored, while at large magnetic fields and chemical poten¬ 
tials the CDW is again favored, this time with large wave 
numbers [43]. 

In summary, the presence of a magnetic field quali¬ 
tatively alters the model phase structure at finite den¬ 
sity, dramatically enhancing the size of the inhomogeneous 
phase. This in turn implies that if quark matter is present 
in the core of compact stars, it might be entirely in a crys¬ 
talline state. Even if for small magnetic fields this would 
only lead to small differences in the equation of state, other 
quantities like transport properties and neutrino emissiv- 
ity might be strongly affected, as we will argue in sect. [^ 

3.3 Vector interactions 

Vector interactions are known to play an important role in 
dense matter. A classical example is the Walecka model jUj, 
where the repulsive coupling of the nucleons to the omega- 
meson field is crucial to reproduce the correct binding en¬ 
ergy of nuclear matter at saturation density. In compact 
stars, where matter is expected to reach several times 
nuclear-matter density, vector-interaction effects should 
therefore be particularly important. 

Vector interactions are known to have a strong in¬ 
fluence on the properties of homogeneous matter in the 
NJL model, where they cause a significant stiffening of the 


Fig. 10. H — H phase diagram allowing for a hybrid chiral 
condensate. Area A corresponds to the CDW phase with q oc 
H, area B is the HCC phase, area C is the large-g CDW phase 
and area D is the chirally restored phase. From ref. m- 

equation of state and reduce the strength of the first-order 
phase transition, possibly resulting in the disappearance 
of the critical point from the phase diagram. 

The effect of vector interactions on inhomogeneous mat¬ 
ter has been investigated in ref. jH]. To this end the NJL 
Lagrangian (eq. 0 ) was extended by the term 

Cv = f , (25) 

where Gv is the vector coupling constant. Its natural 
size, obtained from the Fierz transformation of a color- 
current interaction (“heavy gluon exchange”) is one half 
of the scalar coupling G, but fits to the vacuum mass of 
the omega meson [451H5] or to finite-temperature lattice- 
QCD results at vanishing chemical potential juinsi sug¬ 
gest higher or lower values, respectively. In [2l| Gv was 
therefore treated as a free parameter and varied between 
zero and G. 

Within the mean-field approximation, only the tempo¬ 
ral part V° = (^/’ 7 °' 0 ) of the vector condensate, which is 
identical to the quark number density n, was retainedj^ 
The effect of this field can then be absorbed in a shift of 
the chemical potential, 

/i —>■ /2 = ^ — 2Gvn , (26) 

which is a standard technique for homogeneous matter |5Sj . 
When dealing with inhomogeneous condensates, the quark 
density is in general spatially modulated, implying that 
the shifted chemical potential also depends on the posi¬ 
tion. For a CDW this is not the case, as the density is con¬ 
stant, whereas it is for the RKC modulation (see sect.[^. 

® In homogeneous matter the spatial part V = (ip-fip) of 
the vector condensate does not form, since it would break ro¬ 
tational invariance. This argument does not apply to inhomo¬ 
geneous matter, where rotational symmetry is already broken 
by the scalar and pseudoscalar condensates. Yet, V has also a 
preferred orientation, i. e. it breaks parity, and, at least for the 
RKC, it does not appear |49) . 














































































10 


Michael Buballa, Stefano Carignano: Inhomogeneous chiral symmetry breaking in dense neutron-star matter 




E (MeV/fm^) 


Fig. 11. Equations of state for homogeneous and inhomoge¬ 
neous (CDW and RKC) quark matter in the presence of vector 
interactions with Gv ~ G/2. The inset shows an enlarged de¬ 
tail. 


Nevertheless, in [2T] the shift in the chemical potential 
was implemented only through the spatial average of the 
vector condensate, an approximation which is expected to 
be rather inaccurate close to the onset of the inhomoge¬ 
neous island (cf. fig. [^, but otherwise reasonably valid, 
particularly when approaching the restored phase. 

For Gy = 0, as seen in fig. the inhomogeneous 
phase covers the first-order phase boundary between the 
homogeneous broken and restored phases, and ends at the 
critical point [HUM]. Since vector interactions weaken or 
even remove this homogeneous first-order phase transi¬ 
tion, it was expected that the inhomogeneous phase be¬ 
comes smaller or disappears as well when Gy is increased. 
It turned out, however, that the inhomogeneous phase be¬ 
comes larger instead, keeping its extension in tempera¬ 
ture and enhancing its size in the chemical potential di¬ 
rection m- 


In order to investigate the consequences of the vector 
interaction for the equation of state, we show in fig. [TT] the 
results for homogeneous matter, compared with the RKC 
and CDW modulations, for Gy = G/2. The first visible 
effect compared to fig.j^is the absence of horizontal jumps 
in the curves, as all phase transitions become second-order 
for sufficiently large values of the vector coupling. Addi¬ 
tionally, we observe that at low energy densities the equa¬ 
tions of state for homogeneous and inhomogeneous matter 
cross, i.e. inhomogeneous matter, and in particular the 
RKC, is softer at low and stiffer at high densities. Just 
like for homogeneous matter, the most pronounced effect 
of including vector interactions is in any case the overall 
stiffening of the equation of state, as can be seen in fig. [T^ 
where the pressure with and without vector interactions 
is plotted for the CDW modulation. 

For asymmetric quark matter, the isospin structure of 
the vector interaction becomes again important. In this 
context, instead of the pure isoscalar interaction eq. (25), 
some authors employ the Lagrangian [501IM] 


CyA = -Gy[(V'7^'0)^ + 

+ , (27) 


Fig. 12. Equations of state for quark matter with a CDW 
modulation for Gy = 0 and Gy = G/2. 


which also contains isovector terms and can be viewed 
as t he vector - axial-vector analogue of the Gi term in 
eq. (21). Indeed, as in the latter, the different flavor con¬ 
tributions decouple. Hence, in the mean-field approxima¬ 
tion the effect of the vector part (upper line) reduces to 
a shift of the individual quark chemical potentials for 
each flavor depending on their own number density, i.e. 
Jlj: = — AGyUf. For two-flavor isospin symmetric mat¬ 

ter, where = /id = /i and Uu = rid = n/2, this reduces 
to the standard expression eq. (26). 

Finally we note that, in contrast to eq. (25), CyA con¬ 
tains axial-vector terms. These are required in the isovec¬ 
tor channel by chiral symmetry. Whether the correspond¬ 
ing mean fields have an effect, e.g. on a CDW, has not yet 
been investigated. 


3.4 Competition with color superconductivity 

So far we have concentrated on homogeneous and inhomo¬ 
geneous chiral condensates, i. e. the condensation of quarks 
with antiquarks or quark-holes. It is known however that 
at very high densities particle-hole pairing is disfavored 
against color superconductivity (CSC) [S2], where quarks 
are paired with other quarks. Moreover, NJL-model [551 
I54II55] and Dyson-Schwinger |56j studies of isospin-sym- 
metric homogeneous matter suggest that CSC phases ap¬ 
pear already at the chiral restoration transition. This raises 
the question of how the inhomogeneous chiral condensates, 
which are found in the same region, are affected by CSC, 
and in particular whether they survive at all. 

In refs. [571158] this was investigated within a two-flavor 
NJL model. To this end, the Lagrangian eq. ([^ was aug¬ 
mented by the term 


Cd = Gd (V'^GyV^A^V')^ , (28) 

where X^, A = 2,5,7 denote the antisymmetric Gell-Mann 
matrices in color space and G is the charge conjugation 
matrix. The coupling constant G^ is not strongly con¬ 
strained. Similar to the vector coupling Gy, it is often 
treated as an additional parameter, typically of the order 
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of G or lestQ although slightly larger values have also been 
employed in compact-star phenomenology [5U] . 

Eq. ( p^ describes a quark-quark interaction in the 
spin -0 color-antitriplet channel, which is expected to be 
the most attractive one. The corresponding diquark con¬ 
densate 

(29) 

is an order-parameter for the so-called two-flavor super¬ 
conducting (2SC) phase. It is proportional to the gap pa¬ 
rameter Aj), which can be identified with the gap in the 
quark spectrum if no other condensates are present. 

In refs. [sTiinE] the formation of a nonzero (albeit spa¬ 
tially constant) diquark gap Z\d was studied together with 
a CDW ansatz (eq. (10)) in the chiral sector. Despite this 


complication, the author was still able to find the eigenval¬ 
ues of the effective Dirac Hamiltonian analytically, keeping 
the combined mean-field problem on a manageable level. 
The phase structure was then obtained by minimizing the 
thermodynamic potential with respect to the amplitude 
A and wave number q of the CDW, together with the gap 
parameter Ao- 

For Gd = G12 the results at T = 0 as functions of the 
chemical potential are displayed in fig. One can see 
three different phases, separated by first-order phase tran¬ 
sitions: For fi < 300 MeV, we have Z\ 7 ^ 0 and q = Ao = Q, 
i.e. the normal (non-superconducting) homogeneous chi- 
rally broken phase is favored. For ^ > 400 MeV, on the 
other hand, we have A^ ^ 0 and Z\ = 0, corresponding to 
the 2SC phase with restored chiral symmetry. Finally, the 
most interesting regime is the intermediate one, where the 
CDW (Z\, <7 yf 0) coexists with the CSC gap {Ajj 7 ^ 0). 
Note that there is no region where the CDW exists with¬ 
out CSC. This can be traced back to the fact that atT — 0 
the Cooper instability which leads to Ajj ^ 0 occurs as 
soon as there is a nonzero density. At larger temperatures, 
on the other hand, CDW regions without CSC are possi¬ 
ble [55] . 

From the above argument, it also follows that for small 
values of Gu the location of the coexistence region at 
T = 0 should roughly coincide with the inhomogeneous 
phase in the case without considering CSC|^ Increasing 
Gd favors phases with large diquark gaps. As a conse¬ 
quence all phase boundaries move towards lower chemical 
potentials, albeit in an asymmetric way, so that the size 
of the coexistence phase decreases |57|. 

It is important to note that these studies on the in¬ 
terplay of inhomogeneous chiral symmetry breaking with 
color superconductivity have been performed for isospin 
symmetric matter only, and that the picture described 
here might be significantly altered if this condition is re¬ 
laxed. Indeed, it is known that imposing electric neutral¬ 
ity strongly disfavors 2SC pairing, since the Fermi sur¬ 
faces of the up and down quarks, which are paired in 


^ Performing a Fierz tranformation of a color-current inter¬ 
action yields Gd = 0.75 G. 

® Note that a different regularization was used in ref. [57] , so 
that we cannot directly compare it with the other results we 
have shown. 



0.3 0.35 0.4 0.45 0.5 


Fig. 13. Amplitude A (top) and wave number q (middle) of 
the CDW modulation, and 2SC gap parameter Ad (bottom) 
for Gd = G/2 as functions of p at T = 0. All numbers are 
given in units of GeV. From m- 


the condensate (29), are pulled apart |60|. In particu¬ 
lar, the homogeneous 2SC phase could become unstable 
with respect to the formation of inhomogeneous color¬ 
superconducting phases | 6 T] . Since the corresponding free- 
energy gains are much smaller than for homogeneous BCS 
pairing in isospin symmetric matter, there should be more 
room for inhomogeneous chiral condensates, which, as we 
have seen, are much less sensitive to isospin asymmetry. 


3.5 Strange quarks 

If quarks are present in the cores of compact stars, the 
density might even be high enough to populate strange- 
quark states. This is a rather common scenario in bag 
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models where the strange quark mass is usually taken to 
be of the order of its bare value, < 150 MeV, well below 
the associated chemical potential. Since strange quarks 
are negatively charged, they help neutralizing the matter 
and thus reduce the isospin asymmetry. Moreover, a high 
abundance of strange quarks favors color superconductiv¬ 
ity in the color-flavor locked (CFL) phase where up-, down 
and strange quarks are paired in a symmetric way. This 
phase is electrically neutral even in absence of electrons, 
i.e. for /.tg = 0 [55] ■ 

In the NJL-model, on the other hand, strange quarks 
have considerably higher dynamically masses and there¬ 
fore only appear at relatively high chemical potentials, 
fi > 400 MeV [63] or even 500 MeV [64] Q As a conse¬ 
quence, it is typically found that the densities needed to 
have sizeable amounts of strange quarks (and in particular 
CFL pairing) are either not reached in compact stars, or 
that the softening of the equation of state due to the ap¬ 
pearance of strange quarks immediately renders the star 
unstable [5511551155] . 

Nevertheless, even if their density vanishes, strange- 
quark degrees of freedom can influence the equation of 
state via flavor mixing with the up and down quarks. This 
mixing is caused by the’t Hooft interaction 


Ch = K I det 




det# 




(30) 


which corresponds to the three-flavor analogue to the G 2 - 
term in eq. (21). Here det/ denotes a determinant in flavor 
space, i.e., Ch is a six-point interaction, mixing all three 
flavors. As a consequence the strange quark-condensate, 
which enters the condensate part of the thermodynamic 
potential, is sensitive to the density of the non-strange 
quarks and therefore varies already below the threshold 
for strange quarks. 

Unfortunately, while there exist several publications 
on three-flavor crystalline color superconductors, there are 
so far only very few works dealing with inhomogeneous 
chiral condensates in three-flavor matter. The first and 
most elaborate one is Ref. m, where a three-flavor NJL 
model with ’t Hooft interaction was employed. Within 
that model the authors considered a generalized CDW 
ansatz, keeping the condensate functions eq. (10) in the 
up- and down-quark sector, while the strange-quark con¬ 
densate was taken as a homogeneous background, (ss) = 
const. For this ansatz the mean-field Hamiltonian could 
still be diagonalized analytically, thus keeping the prob¬ 
lem on a tractable level. Nevertheless, since the strange 
and non-strange sectors are coupled through the’t Hooft 
term, the strange-quark condensate has a non-trivial ef¬ 
fect on the phase structure. In particular it was found that 
the chemical-potential window where the inhomogeneous 
solution is favored can increase considerably. 

The calculations of Ref. m have been performed for 
a common quark chemical potential. It is likely that a 


® Similar thresholds have also been found in Dyson- 
Schwinger QCD studies of color-superconducting matter with¬ 
out neutrality constraints [56) . 


similar enhancement of the inhomogeneous phase will also 
occur when neutrality conditions are imposed, although 
this still needs to be checked explicitly. 

In Ref. [68] the same CDW-like modulation was con¬ 
sidered in electrically neutral matter but without’t Hooft 
interaction. The influence of the strange quarks on the 
non-strange sector in this case is therefore only indirect 
via the neutrality condition, and hence, their effect on the 
inhomogeneous phase is quite marginal, unless very high 
densities are reached in the core of compact stars. 


4 Astrophysical implications of 
inhomogeneous chiral symmetry breaking 


After describing the thermodynamical properties of inho¬ 
mogeneous chiral symmetry breaking phases within effec¬ 
tive models of cold and dense quark matter, we now want 
to discuss possible consequences of the formation of these 
crystalline structures on phenomenological observables re¬ 
lated to compact stars. 

The determination of clear signatures relating macro¬ 
scopic observables with the microscopic properties of mat¬ 
ter within stellar objects is obviously a very challenging 
task. From the results for the equations of state we have 
discussed earlier we can already anticipate that the effect 
of inhomogeneous condensates on the mass-radius relation 
will be marginal. This will be discussed in sect. |4.l| More 
promising signatures can be expected from observables re¬ 
lated to transport properties, since they are directly in¬ 
fluenced by the crystalline structure of the matter and 
the corresponding low-energy excitations (phonons). Some 
of these signatures, like gravitational waves and glitches, 
have been investigated for crystalline color superconduc¬ 
tors (cf. ref. [H]), but not yet for the specific case of in¬ 
homogeneous chiral condensates. An exception is the neu¬ 
trino emissivity, which has been worked out for a CDW in 
ref. [59] • This will be reviewed in sect. 4.2 


4.1 Mass-radius relations 

A possible way of discerning between different microscopic 
compositions of matter in a compact star is to look at 
the relationship between its mass and radius. Mass-radius 
sequences relating these two quantities can be obtained 
by solving the Tolman-Oppenheimer-Volkoff (TOV) equa¬ 
tion [Toiin], which uses the equation of state of the sys¬ 
tem as an input. Until relatively recently, this method has 
been unable to provide clear information, as simultaneous 
measurements of masses and radii of compact stars are 
extremely challenging and plagued with large uncertain¬ 
ties. The situation partially changed in the past few years 
thanks to the very precise measurements of pulsars with 
unusually large masses, close to 2Mq mil], Mq being the 
solar mass. These measurements allowed to rule out a large 
number of softer equations of state, which would be unable 
to support such massive stars, and serve as a benchmark 
for all possible models describing the microscopic compo¬ 
sition of matter inside these compact objects. 
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Fig. 14. Mass-radius sequences for a pure quark star, neglect¬ 
ing electric neutrality and magnetic fields. Top: Comparison of 
three different modulations for Gv ~ 0. Bottom: CDW modu¬ 
lation with different strengths of the vector coupling. 


As discussed in sects. and the appearance of in¬ 
homogeneous condensates results in relatively small vari¬ 
ations of the equation of state, especially when compared 
with the uncertainties caused by using different models 
and parametrizations. As a consequence, the effects of 
spatial modulations on the mass-radius relations are very 
small as well. To see this, we show in fig. [^mass-radius se¬ 
quences of different equations of state. Since at this point, 
rather than providing a realistic description of stellar mat¬ 
ter, we just want to compare the influence of different 
mechanisms, we consider pure quark stars and neglect 
electric neutrality and magnetic fields. In the upper panel 
the equations of state of fig.[^have been used, correspond¬ 
ing to different spatial modulations in matter without vec¬ 
tor interactions!^ Except for the low-mass region, which 
should become irrelevant once a more realistic description 
of the hadronic phase is implemented, the effect of the 
inhomogeneity is negligible. Comparing this result with 
the curves shown in the lower panel, which have all been 
obtained with the CDW ansatz, but for different values 
of the vector coupling, we conclude that it would be im¬ 
possible to infer the existence of an inhomogeneous phase 


Note that these equations of state do not correspond to 
self-bound quark matter. Therefore the curves turn around to¬ 
wards large radii at low masses, similar to ordinary neutron 
stars, but unlike strange stars in the case of absolutely stable 
strange quark matter. 


on the basis of a mass-radius measurement, even if it was 
very precise. 

On the other hand we have seen in sect. El that the in¬ 
homogeneous phase itself is rather robust and that its size 
is even increased by the presence of vector interactions, 
strange-quark degrees of freedom, and magnetic fields. 
This means that, as long as there is a quark phase at all in 
compact stars, it is quite likely that it is inhomogeneous [^ 

In order to determine whether a compact star with 
a crystalline core is compatible with the 2Mq measure¬ 
ments, the authors of [55] investigated solutions of the 
TOV equation for a hybrid configuration, consisting of 
an inhomogeneous quark-matter core surrounded by nu¬ 
clear matter. The latter was described using the non-linear 
Walecka model, while for the quark phase the authors 
considered a CDW modulation within an extended NJL 
model. Vector interactions, electric neutrality and beta 
equilibrium were implemented, considering for simplicity 
a Lagrangian without flavor mixing, i.e. eq. (21) with 
G 2 = 0 for the scalar-pseudoscalar interaction and^q. (27) 
for the vector interaction. The effect of strangeness was 
also studied, employing three-flavor versions of the NJL 
model (without’t Hooft interaction) and of the non-linear 
Walecka model. Following ref. 162] inhomogeneous quark 
matter was then set up as a CDW in the up- and down 
quark sectors in the homogeneous background of a spa¬ 
tially constant strange-quark condensate. Color supercon¬ 
ductivity, on the other hand, was not considered. Finally, a 
magnetic field with a chemical-potential dependent strength 
varying from lO^^G at the surface to 2.5 x 10^®G in the 
core was included in the model. 

Within this approach, hybrid stars with maximum mass¬ 
es above 2 Mq and containing a quark-matter core in the 
inhomogeneous phase have been obtained, suggesting that 
such configurations are possible in nature. It is important 
however to recall that in the presence of magnetic fields, 
neutron stars are no longer spherical. Therefore the use of 
the TOV equations in this case is incorrect, and a calcula¬ 
tion which properly takes into account the deformed star 
shape is required Nevertheless, the expected 

influence of the magnetic field on the maximum mass is 
relatively small [33], so that a more refined calculation 
would likely lead to values which are compatible with the 
ones of [55] , especially if one considers the large uncertain¬ 
ties originating from the use of phenomenological models 
with free parameters. A numerical cross-check of this ex¬ 
pectation would be of course extremely desirable. 

A major uncertainty in modeling hybrid stars arises 
from the fact that usually two different models are em¬ 
ployed to describe hadron and quark phases. As a conse¬ 
quence there is in principle an unknown parameter (“bag 
constant”) which determines the relative difference of the 
pressure in the two models m- In addition, there are dif¬ 
ferent ways how to join them. While in most calculations, 
including Ref. [68], a Maxwell construction is performed, 
it was argued in Ref. m that the more correct approach 


If high enough densities are reached without having a 
phase transition to quark matter, inhomogeneous chiral con¬ 
densates could also become favored in the hadronic phase m- 
















14 


Michael Buballa, Stefano Carignano: Inhomogeneous chiral symmetry breaking in dense neutron-star matter 


is a Gibbs construction. This is true, however, only for 
sufficiently small values of the surface tension, which is 
practically unknown, thus being another source of uncer¬ 
tainties. 

Ultimately, one should therefore aim for a unified de¬ 
scription of the quark and hadron phases. A model allow¬ 
ing for inhomogeneous quark matter might be a step in the 
right direction [32], as it could allow for quarks to lump 
together into baryons. Indeed, it is known since more than 
25 years that baryons in the NJL model can be described 
as non-topological solitons [THirffllTH] . which are associ¬ 
ated with a spherically symmetric inhomogeneous chiral 
field. Obviously it would be very exciting to pursue this 
idea further. 


chiral restoration transition, where the deformation is in¬ 
stead very pronounced. In the first case, the result for the 
emissivity is given by 


onset 

^CDW 
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while in the second it is 
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4.2 Neutrino emissivity 

A perhaps more promising way to get insights about the 
composition of compact stars is to observe their cooling 
process, which is expected to happen mainly through neu¬ 
trino emission. A correct phenomenological description of 
this process should be able to account for the relatively 
low surface temperatures of young pulsars like the Vela or 
3C58 [751150] . without conflicting with the slower cooling 
of Gas A [5T|. In particular, since the mass of Gas A was 
estimated to be > 1.5M0, the cooling rate should not rise 
too fast with the mass of the star. This suggests that some 
exotic cooling scenario might be required to explain these 
results. 

Among these mechanisms, quark cooling is very ef¬ 
ficient, but the direct Urea processes d ^ u + e~ + Ve 
and u + e~ —>■ d are kinematically suppressed for 
non-interacting massless quarks. This is encoded in the 
so-called triangular condition, a direct consequence of en¬ 
ergy and momentum conservation stating that in this case 
all momenta must be collinear, thus dramatically reduc¬ 
ing the phase space available. This in turn implies that 
in the standard scenario of a high-density phase of chi- 
rally restored quarks, at least for low temperatures, these 
processes cannot account for the very rapid cooling ob¬ 
served. The picture however changes significantly when 
inhomogeneous phases are included: if a crystalline chi¬ 
ral condensate forms, quark momenta will no longer be 
conserved due to their interactions with it. The inhomo¬ 
geneous condensate therefore provides the additional mo¬ 
mentum which might relax the triangular condition and 
make these processes kinematically allowed, without the 
need for spectator particles. 

This scenario has been worked out in detail in [55] , 
where the neutrino emissivity (ie. the emitted energy per 
unit volume and time) has been calculated for the case of 
a GDW modulation (eq. @) within a simple two-flavor 
NJL model. In the following we will report the main re¬ 
sults of their calculation. Since a full computation becomes 
extremely involved due to the deformation of the Fermi 
surface generated by the formation of inhomogeneous con¬ 
densates, the authors focused on two limiting cases: the 
onset of the inhomogeneous phase, where the Fermi sphere 
deformation is relatively small, and the region close to the 


Here Gp = Gp cos 9c, with Gp weak coupling constant 
and 9c Cabibbo angle. For the first case, the dependence 
on the wave number of the GDW modulation is of the 
form 1/g, arising from the integration over the phase fac¬ 
tor associated with the inhomogeneous chiral con¬ 

densate. This specific momentum dependence is similar 
to the one obtained for pion condensation, and the re¬ 
sulting magnitude of the emissivity is comparable to the 
efficient one obtained for massive quarks. Glose to the 
restoration transition the g-dependence becomes instead 
more involved, and the emissivity is typically up to two 
orders of magnitude lower. At any rate, both results are 
significantly larger than standard processes like the mod¬ 
ified Urea. This allows for quick cooling, similar to the 
Urea process with interacting quarks [53] or other exotic 
mechanisms such as the ones related to the presence of a 
pion condensate [55] . 

While this calculation has been performed only for the 
specific case of a GDW modulation, we can expect neu¬ 
trino emission to be significantly enhanced for any kind of 
modulation of the chiral condensate. A similar calculation 
for the RKG would be desirable, but likely to be more in¬ 
volved due to the more complicated structure of the quark 
eigenfunctions required to calculate the transition matrix 
elements. 

Ultimately, the exact cooling rate of a star will of 
course depend on the extension of the inhomogeneous win¬ 
dow. This in turn implies that cooling measurements might 
in principle impose limitations on the size of inhomoge¬ 
neous phases: if crystalline chiral condensates form only 
in a limited range of densities, cooling will be most effi¬ 
cient within a relatively thin shell inside the star. Hence, 
as long as chiral restoration is reached in the center of 
the star, the emitted energy grows only like the surface of 
that shell. The authors of ref. [55] therefore conclude that 
heavier stars may not cool faster than lighter ones. 

Gonversely, if the inhomogeneous phase is very broad, 
more massive stars will have a larger region occupied by 
crystalline condensates and will cool more efficiently. In 
order to support these qualitative arguments with a more 
quantitative estimate, however, more thorough investiga¬ 
tions are needed, most notably an evaluation of the heat 
capacity of the inhomogeneous layer, as well as its thermal 
coupling to the other regions of the star. 
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5 Conclusions 

We discussed the phenomenon of inhomogeneous chiral 
symmetry breaking, with a particular focus on its astro- 
physical implications. Several models suggest that the for¬ 
mation of spatially modulated structures of the chiral or¬ 
der parameters is energetically favored in a window of the 
phase diagram which is relevant for cold and dense stel¬ 
lar matter. After introducing some basic thermodynami¬ 
cal properties of these crystalline solutions, we discussed 
how they are affected by including several necessary con¬ 
straints for a correct description of matter in compact 
stars. While electric neutrality and beta equilibrium do 
not strongly affect the size of the inhomogeneous window, 
other effects, like vector interactions or the presence of 
background magnetic fields, enlarge it significantly. The 
role of color superconductivity, on the other hand, is not 
yet completely clear. According to NJL-model studies in 
isospin symmetric matter inhomogeneous chiral con¬ 
densates may coexist with a homogeneous diquark con¬ 
densate if the diquark coupling is small, but become sup¬ 
pressed for larger couplings. This could change, however, 
if electric neutrality is imposed, since it strongly disfavors 
the BCS pairing of up and down quarks. As a possible re¬ 
sult, the diquark condensates could become spatially mod¬ 
ulated as well or be suppressed completely. In any case, 
this implies that, if quark matter is present inside the core 
of these stars, it is most likely in a crystalline state. In par¬ 
ticular, the realization of a coexistence phase where both 
color-superconducting and chiral condensates are inhomo¬ 
geneous is a fascinating possibility, which might be worth 
investigating. 

Hybrid stars with an inhomogeneous quark matter- 
core are compatible (within model uncertainties) with the 
recent measurements of pulsars with 2 solar masses |68] . 
Since the presence of crystalline phases has however only 
a small influence on the equation of state, their existence 
or non-existence cannot be concluded on the basis of mass 
measurements, even in conjunction with precisely deter¬ 
mined compact-star radii. Indeed, model uncertainties and 
systematic errors such as those introduced by employing 
spherically-symmetric TOV equations in an anisotropic 
system are expected to dominate over these effects within 
this type of calculation. It is therefore highly desirable to 
work out more sensitive signatures of their formation. A 
first example in this direction is the neutrino emissivity, 
discussed in , which is related to the cooling properties 
of the star. A complete calculation of the cooling process 
of the whole star would however be required before a de¬ 
tailed comparison with astrophysical data could be made. 

The formation of inhomogeneous phases might also 
have a stronger influence on the transport properties of 
matter inside compact stars, particularly through their pe¬ 
culiar low-energy excitations, such as phonons. Some work 
on calculating the dispersion relations for these excitations 
has been performed in [84] for the case of a CDW mod¬ 
ulation and [SS] for a RKC. It would also be interesting 
to investigate phenomena like glitches and the emission of 
gravitational waves, which have been investigated for crys¬ 


talline color superconductors (cf. ref. (Hj), but not yet for 
the specific case of inhomogeneous chiral condensates. 
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